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ABSTRACT 

Most physical systems are nonlinear to some extent. One variety 
of nonlinear system uses a relay or other on-off device as a control 
element. 

Since in the process industries or relay servo systems, one often 
encounters a type of time delay called transportation lag or dead time. 

It is highly desirable to investigate the effects of the transportation 
lag on the process or relay servo system. 

This paper is to demonstrate the effects of the transportation lag 
on the relay servo system and show how to eliminate these effects. 

Experimental results indicate that the transportation lag causes a 
continuous oscillation with a steady state over- and under-shoot, the 
response can never be settled, and also it leads to a longer rise time. 

In addition to those, the transportation lag results in a terminal limit 
cycle, Though the tachometer feedback has been used to compensate the 
system, this oscillation can never be eliminated completely when some delay 
exists. Therefore, a five percent overshoot and undershoot was considered 
acceptable. 

The author wishes to express his appreciation for the advice, assist- 
ance and encouragement given by Dr. Richard Carl Dorf, Professor of Electri- 
cal Engineering, United States Naval Postgraduate School, throughout this 
mest ization; He also owes a debt of gratitude to Dr. John R. Ward, 
Professor of Electrical Engineering Department, U. S. Naval Postgraduate 


School, for his assistance during the initial stages of the work. 
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1. Introduction, 

Systems which utilize feedback for control purpose have become es- 
sential elements in modern technology. They range all the way from 
simple toys to the most complex automatic factories and production equip- 
ment. Feedback control is the guidance that makes modern automation pos- 
sible, and as such is in a large measure responsible for the ever-increas- 
ing productivity and rising standard of living of man. 

It is noted that feedback control systems are seldom truly linear, 
but rather they are subject to the introduction of a considerable number 
of nonlinearities. Some nonlinearities may be introduced intentionally by 
the system designer, An example of this might be the selection of a re- 
lay servomechanism with its attendent nonlinearity. On the other hand 
there may be no way of avoiding a nonlinearity in the system design. An 
example of this would be the mandatory inclusion of a gear train with its 
attendent backlash. A nonlinearity which frequently can neither be avoided 
nor neglected is transportation lag. 

In some control systems, especially in the field of process control, 
a non-eligible propagation time is encountered, That is, there is a definite 
lapse of time between the initiation of a change in the signal at one point 
in the system and the appearance at some other point of any effect of this 
change. But be sure that this kind of lag is not to be confused with the 
type of lag introduced by factors of the form 1/(TS + 1) or Wits? + 2 3° TS 
+1) in the transfer function relating the two signals. Factors such as these 
in the transfer function of a component do reveal a tendency of the com- 
ponent to retard the signal, in a sense, by introducing phase lag at all 
frequencies, but they do not result in a finite dead time which must elapse 
before any effect of a change in input can be noted in the output of the 


component, 





Lead time of appreciable magnitude can be observed, for example in 
long pipes carrying fluid under pressure. A change in pressure at one end 
of such a pipe is transmitted with finite velocity to the other end, pos- 
sibly but not necessarily being distorted in transmission. Another type 
of system in which dead time is often encountered is shown in Fig. 1-1. 

In such a system two different fluids are mixed at the mixing point, and 
the relative concentration is read at the metering point. If the meter- 
ing point is located at an appreciable distance from the mixing point, 
the dead time between the establishing of a change in relative concentra- 
tion and the detection of the change may be appreciable. 

Another example is illustrated below. In Fig. 1-2 the shaded area 
represents a slug of water at position 1, the outlet of the steam nozzle. 
Position 2 represents where the temperature sensing element will be placed. 
Let the temperature of the shaded area water slug at the steam nozzle be 
quickly changed from 6, to Oe . Before the change can be measured, the 
slug of water must travel from the steam nozzle to position 2 where the 
sensing element would be located. 

The effect of the delay caused by the slug traveling down the pipe 
can be shown by locating sensing elements at position 1 and 2 as shown in 
Fig. 1-3. If the time constants of the sensing element at position 1 is 
taken as the input and the element at position 2 as the output, Fig. 1-4 
becomes the reaction curve for the process. 

This type of reaction curve is typical of a dead time lag (or trans- 
portation lag) and may be defined by the time lag T. The magnitude of 
the dead time lag will depend on the water velocity and the distance be- 


tween position 1 and 2, The lag time may be calculated as follows: 


T(min.) = d(ft) / V(C£t/min) 
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where, T = dead time lag (min) 


d = length (ft) ~ 


" 


V = velocity of flowing medium (ft/min) 

Dead time lag is very common in the process industries. Because 
of the complete lack of the response during the dead time, this type of 
lag is very detrimental to control. An examination of the equation de- 
fining the dead time lag shows that increasing the length or decreasing the 
velocity will increase the lag time, and we shall see later, make the 
problem of control more difficult. 

To determine the frequency response of the dead time lag, a sinuscilal 
wave having an amplitude of 10°F will be generated around a bias tempera- 
ture o1° At some low frequency the curve as recorded by the sensing ele- 
ments at position 1 and 2 would appear as shown in Fig. 1-5(a). 

Fig. 1-5(b) and 5({c) show similar results for two successively higher 
frequencies, As can be seen, as the frequency increases the period decreas- 
es. However, the lag time which is dependent only on the velocity and 
length remains constant. Thus the magnitude of the phase angle will in- 
crease proportionally to each decrease in period (increase in frequency). 


The phase angle may be calculated as follows: 


§ (degree) eee ts) X 360 (degree/cycle) 


Period (min/cycle) 
or ® (degree) = - T(min) X f(cpm) XK 360 (degree/cycle) 


From this equation we see that the phase angle for the dead time lag 
is directly proportional to the frequency. The phase angle approaches 
zero degrees as the frequency approaches zero and negative infinity as 


the frequency approaches infinity. This will be discussed in © 6. 
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2. Relay servo system, 


2.1 General Description. 

When economy is of importance, many feedback control systems utilize 
a contactor for energizing the actuator. Such a system is called a contac- 
tor or, alternatively, an on-off control system, since the contactor is an 
element whose output can usually have only two values. In some cases 
these correspond to a fixed positive action and no actuator action. This 
usually is the case with thermostatic controls, for example. In other 
cases the two values correspond to maximum positive and negative actuator 
action. 

In Fig. 2-1 it is assumed that N is a relay or switch which has an 
output +1 or -1, depending on the sign of E = R- b. A positive output cor- 
responds to a positive sign of E and a negative output corresponds to a 
negative sign of E. When the switch is a perfect switch this can be de- 
noted symbolically by the equation 


V = Sign E (uals) 





Fig. 2-1 Block diagram cf the contactor servomechanism 
Since the contactor is frequently an electrically controlled or equivalent 
element, there may be a delay in the initiation of switching action as 
well as a finite time interval during which switching action takes place. 
The contactor action ordinarily applies full control force or torque to 


the system. Since the output of N is therefore not necessarily unity, the 





actual value is included in the gain constant of the following element. 
The block G in Fig. 2-1 will be assumed to be a linear servomotor 
with viscous damping and no load torque except that due to inertia. The 
limiting case of zero viscous damping is of interest and will also be 
considered, Including the gain constant of the switch, the motor transfer 
function is 
fe 
FSCS * Ven) 
where To = KE | is the stall torque of the motor and J is the total inertia 


G = (252) 
of motor and load. In the absence of damping, l/r is set to zero in this 
equation. It should be emphasized that since N is nonlinear element, there 
is no particular reason to assume that the feedback element is linear. [In 
particular, the feedback signal could be any algebraic function of position 
and its derivatives, However, for the method of analysis used here the 
feedback function can depend only on the position and velocity and must 
not on the derivative of the velocity (acceleration). The symbol H, there- 
fore, no longer represents a linear transfer function, but a possibility 
nonlinear feedback function involving the output position and velocity. 
Symbolically this will be written as 

b = £(C,C) (2,3) 
where C = dC/dt 
2.2 Relay servos 


The system behavior can be prescribed by combining the equations: 





ve i 2 Aceren To >: is ~ 
oe, oO FT eee) Se 
“ef 3 o vr) e > : e ie 
x tHe =m— —L2 as (c)=< © 
Zo» JS 


Between switching times the response is described by a second order 


differential equation. For this reason such a system is called a second- 





order contactor servo or feedback system. The two equations are inte- 
grable between switching times and the entire solution can be obtained 
by matching final and initial conditions of the two sets of solutions at 
the proper switching time. The solution is therefore easy to obtain in 
principle. The details, however, are so complex that another method of 
solution is desirable. 

The analysis of the system can be simplified by using the fact that 
it is possible to plot the response of the system as a trajectory in the 
phase plane, the trajectory being composed of a series of arcs or curves 
which have a characteristic shape dependent upon the properties of the 
system, The behavior of a physical system where only one coordinate is of 
interest is usually specified by determining this coordinate as a function 
of time, C = fy (t). This function depends upon the differential equation 
describing the system and upon a set of initial conditions. For some 
purpose, an equal suitable way of describing the behavior of the system 


is by means of a functional relationship between position and velocity: 


C=__ f, (c) C25) 
When such a4 relation is known, the time dependence of C, if required, 


can be found directly by integrating the equation 





re a, 
: fae: ae | . SPL 
a /, Co wes ae =f Acc) (2.6) 


Equation 2.5 can be represented graphically with position and velocity 
taken as rectangular cartesian coordinates; position is usually plotted 
as abscissa and velocity as ordinate. The behavior of the system is des- 
cribed by a plot of all successive pairs of values of position and velo- 


city. The plane on which this plot is made is called the phase plane, and 





the graph of Equation 2,5 is called the phase-plane trajectory. 

For the system to be considered in this paper, the differential 
equations of interest will not contain time explicitly. For the 
second-order system this implies that the acceleration can be written as 

C =F (ce, cy 
The importance of this concept here results from the fact that for such 
systems, the behavior and therefore the trajectory are uniquely specified by 
the initial values of C and C. This means that there is an only tra- 
jectory through a given (nonsingular) point in the phase plane. [In this 
paper it is restricted to the systems described by second order differential 
equations and hence by phase plane trajectory in the two dimensional C= 
phase-plane. 


The slope of the trajectory in the phase plane can be written as 


Jc 
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According to this equation the slope is equal to the quotient of ac- 
celeration and velocity. In the upper half plane, C is positive and 
nonzero, and therefore the slope is nonifinite and either positive nega- 
tive, depending upon whether the acceleration is positive or negative. 
Similarly, in the lower half plane the slope is non-infinite and either 
positive or negative, depending on whether the accele)ation is negative 
or positive. For a point on the horizontal axis the velocity is zero, 
and therefore unless the accelerations also zero the slope dC/dc is in- 
finite. In general, then, trajectories cross the horizontal axis verti- 
cally. Direction of motion along the trajectory can always be determined 


by the fact that for C positive C must become more positive with time and 





for C negative C must become more negative with time. Motion along 
the trajectory is therefore to the right in the upper-half plane and to 
the left in the lower-half plane. The horizontal component of the 
trajectory velocity is equal to the system velocity C = éC/dt and 
the vertical component of the trajectory velocity is equal to the system 
acceleration C = dc/dt. These considerations make it clear that a trajec- 
tory can not end except on the horizontal axis. If there are points on 
the horizontal axis where both C and C are zero they are called 
singular points, and the value of C at such a point corresponds to an 
equilibrium condition of the system. Trajectories can emerge from or 
return to some of these singular points with a non-infinite slope, since 
the slope dC/dC is undefined at this point except through other considera- 
Erons . 
2.3 Phase plane trajectory. 

Some of the ideas in 52,2 will now be illustrated with the system 
shown in Fig. 2.1 and defined by Equation 2.2 through 2.4. The Equation 


2.4, which are second-order differential equations, can be reduced to 


pairs if first-order equations by substitution of Equation 2.6 
r= C = 5G 
— WL 
which gives 


IV / re 18 
het ienaan a ed rea £ = ——- a ? m 


an alternative expression for velocity, dV/dt = V(dV/dC), can be used 


with Equation 2.8 to give the following equation: 


7 fv VY ye 
y+ = - _ Raa (2.9) 


This equation can now be integrated with the result 


——= y, | ~ bile C // SV. 
C-—-C. = -Z,[V 6: 2 lo3(/ ie a =) (2,10) 
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Evidently the constant of the integration Co is the value of C at 
which V = 0. This equation can be put in dimensionless form by the 


substitutions 


- i = y = 
y=-—/, C == -= Cc (2.11) 
ae bin es 





a 


which gives 
C-2, = - [Ut Ce (FV) ] (2.12) 


Typical trajectories for particular values of C. are shown in Fig. 2.2. 
Through use of the dimensionless variables v and c, all systems of this 
type have identical trajectories. 

In the absence of damping, the trajectories are parabolas, If 1/z_, 


is set to zero, integration of Equation 2.8 gives the result 


| ee eel “i, Ve 
Cen —- = a (2,13) 


In this case the substitution 2.11 can not be used. For the case of 
finite damping it is sometimes convenient to have an expression for the 
trajectory which shows the deviation from the parabolic shape more 
directly than Equation 2.10 or 2,12. This can be obtained by use of the 


power series expansion for log (1 + x): 
Bed: 


= i A We = JV = JV et 
cies + a hea TZ (f5) 2 J ele) 
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or 


ee ee Vi Hema) 2 75) 
C C, =i> [iz = V Ft a (V)+ ~=(Y) foes | (2.15) 


These equations are particularly useful when an approximate result is 
desired, For accurate computations Equation 2.10 or 2.12 is more use- 
ful, since the series given above converges rather slowly except for 
small values of v, and indeed diverges for lvl > ja 


Reference to Fig. 2.1 and Equations 2.1 and 2.3 shows that the 


iat 








Fig. 2.2. Phase plane trajectory for the second-order 
contactor servomechanism. 


condition governing switching from one trajectory family to the other is 
given by the equation 

a (Ce ae (2,16) 
For the present it will be convenient to restrict attention to step in- 
puts with initial velocity zero. In this case a simple translation of 


the origin can be used so that the initial conditions are c 


i 


“Cys Vea. 
and the final desired steady-state conditions are R=C=V= 0. Setting 


R = 0 in Equation 2.16 gives the equation 


#(€,¢)==0O Cay 
The graphical plot of this equation separates the phase plane into two 
parts, in one of which the system follows trajectories corresponding to 
V = +1, and in the other of which the system follow trajectories corres- 
ponding to V = -l, this is illustrated in Fig. 2.3. 
For certain functions £(c,c) an anomalous situation arises. Consider, 
for example, the point P, in Fig. 2.3, which will supposed to be reached 


1 


by motion along the trajectory marked (1). When the switching line is 
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Fig. 2-3. Switching line for contactor servomechanism. 
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reached, the system will normally transfer to the trajectory of the 


other family which pass through P,, this is the trajectory (2) in Fig. 


12 


2.3. However, when the point P, is reached, for example along trajectory 


Z 
(3), a paradoxal situation arises, Although switching should take place 
at this time, there is no trajectory of the other family along which the 
system can move, since P, is an end point rather than a starting point for 
trajectory (4). In terms of real physical systems the paradox can be 
resolved by recognizing that a small delaying switching is unavoidable. 
This permits the system to overshoot the switching line so that a small 
section of allowable trajectory is available in other family. When this 
trajectory reaches the switching line another slight delay allows over- 
shoot back into the left-half plane. The resultant trajectory consists 
of a series of short arcs in which the system oscillates about the switch- 
ing line as it approaches the equilibrium point at the origin. In terms 
of the ideal system, one can say that the trajectory follows the switching 
line to the origin accomplished by an infinite number of switching which 
provide, on the average, the proper force to give an acceleration consist- 
ent with the location in the phase plane and the slope of the switching 
line. 

Using the above considerations some interesting general conclusions 
can be drawn: 

(1) For no viscous damping ( 2, = ©<) and no rate feedback 

[fee cy=2c] the trajectories are parabolas separated by the vertical line 
¢ = 0. The system is neutrally stable in the sense that for any amplitude 
of step input the oscillation will persist indefinitely, as shown in Fig. 


Zee . 


(2) For viscous damping and no rate feedback Fe e)=e |, the 
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(c) Slope = - nai 





trajectory are the “asymmetric parabolas" of Fig. 2.2, separated by the 
verticle line c = 0. The system exhibits a damping or decrement of the 
oscillations produced by a step input. The per-unit decrement is a func- 
tion of amplitude and approaches zero as the amplitude approaches zero, 
In other words, for large input the initial per-unit decrement is high, but 
decreases to zero as the amplitude of the oscillation approaches zero. 
This can be observed qualitatively in Fig. 2.4b. 
(3) With linear rate feedback, the feedback function is fe <)=C*RV 
If there is no viscous damping ( Z,, = 20) , then the trajectories are 
parabolas separated by the straight line 
: / 

— Ke # (2.18) 
Again the system shows a per-unit decrement or damping which varies with 
amplitude. However, in this case the damping is small for large amplitude 
and as the amplitude decreases the damping increases until the final tra- 
jectory into the origin lies along the line 


{ 
Ke 


with no further overshoots. An illustration of this is shown in 2.4c. 
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The above conclusions shows that no simple linear feedback function 
will give optimum response, regardiess of the system's viscous damping. 
If rate feedback is adjusted to give no overshoot for the largest expected 
input then the system moves to the origin slowly for small inputs. If the 
rate feedback is adjusted to give faster response for small input, then 
overshoot and oscillations will follow application of a large step input. 

It is natural to ask whether any feedback function can be found to 
give optimum response for all inputs. This turns out to be a difficult 
question. Since the system is nonlinear and does not have an easily deter- 


mined response for arbitary inputs. However, a definite answer can be found 
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for the case of step inputs by following chain of reasoning. For most 
rapid response without overshoot it is necessary to apply full positive 
acceleration torque until that instant at which full negative torque will 
bring the system to rest at the origin in the phase plane. This is il- 
lustrated in Fig. 2.5 for a particular initial starting point. The tra- 
jectory shown passing through the origin can be designated as the homing 
trajectory. In order to achieve switching on this trajectory for all in- 
puts, it is necessary that feedback function set equal to zero, 

fF Ce. c)=—O0O (2.19) 
being the equation of this homing trajectory. In this case, when any 
trajectory crosses this line switching will occur, so that the applica- 
tion of full reverse torque will just bring the system to reset at the 


origin. 


“(e) 


Homing trajectory 





Fig. 2.5 Switching point for optimum step response 
Since the trajectory of Equation 2.8 is roughly parabolic, a suit- 


able form of the feedback function is 





# C<,v) 2S 1 Se ges 7\¥|)— C+ a y |v (2,20) 
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Strictly speaking, this equation should be modified to match the appropri- 
ate “‘asmmetric parabola'', although for practical purposes the refinement is 
not necessary. In other words, if the position feedback is linear and the 
rate feedback roughly quadratic (except for sign) then it is possible to 
obtain optimum response for this particular type of input. It will be 

seen that switching delays give rise to oscillatory motion of small ampli- 
tude unless a linear velocity-feedback term is introduced. For this reason 
it is usually desirable to introduce a small amount of such feedback (in 
addition to the quadratic feedback) in order to eliminate or minimize this 
oscillation. This does not significantly affect the results for large 
amplitudes. 

It is quite important to remember that this optimization does not 
necessarily give optimum response for arbitary inputs. Further considera- 
tion is beyond the intended scope of this paper. 

2.4 On-off control. 

There probably are more on-off discontinuous controllers in the world 
than the number of all other servomechanisms and regulators put together. 
Hot-water heaters, ovens, refrigerators, furnaces, electric irons, battery 
chargers, liquid-level controls, and compressed-air pressure regulators 
are all commonly of this type. 

A relay system is one which contains a decision function. [t is 
commonly in the error amplifier and in its most elementary form is an on- 
off switch. If the positive error exceeds a certain amount, the switch is 
turned from on to off, and if the negative error exceeds a certain amount, 
the swtich is turned from off to on. When the switch is on, a predetermined 
constant corrective force is applied into the rest of the system to bring 


the output back into correspondence with the desired value. When the switch 


18 





is off, this force is removed. The differential equation for each is the 
open-loop equation, with one term, the corrective force, having two differ- 
ent values for the two periods. Each switch operation starts a new period 
for which the initial values are the final values of the previous period. 

One characteristic of a relay is that it delivers an output having one 
of only two values, independent of the continuous changes in the input, 
except they specify the switching instants, It also has hysteresis, the 
phenomenon that the input values for switching on and for switching off 
are not exactly equal, The actual time of motion of the armature will be 
considered negligible, but the system may contain the dead time of flow of 
material through a long pipe, equivalent to the time of transmission of an 
electric signal on a long distoritonless line. 

The term relay is used to apply to any device having at least the first 
of these characteristics, In this paper, we are interested primarily in 
discontinuous control of continuous processes, so that the relay is only 
one part of the system, 

The on-off control, or two position control, is undoubtedly the most 
widely used type of control for both industrial and domestic service. It 
is the kind of control generally employed on home heating systems and 
domestic water heaters. 

Two position control is a position type of controller action in which 
the manupulated variable is quickly changed to either a maximum minimum 
value depending upon whether the controlled variable is greater or less 
than the set point. The minimum value of the manupulated variable is usually 
zero (off). This mode of control is illustrated by the electric level- 
control in Fig. 2.6. A float in the vessel operates an electric switch 


which controls power to a solenoid valve. When the liquid level rises, 


1s, 








Fig. 2-6, Two position control 
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the switch contacts are closed, the solenoid valve closes, and the flow is 
cut off. When the liquid level falls, the switch contacts are opened, the 
solenoid valve opens, and the inflow resumes. If the float lever has no 
bearing friction and the electrical contact draw no arc, the action is sharp 
or 'knife'edge" as shown in Fig. 2-7. 


The equation for the two-position control are 


m = My when e => 0 
n= M, when e <~ 0 
where m = manupulated variable 
My = maximum value of manupulated variable (on) 
Ea minimum value of manupulated variable (off) 


e = deviation 
Thus two-position control must be described by two equations, each apply- 
ing in a certain region of deviation. 

A differential gap in two-position control causes the manupulated vari- 
able to maintain its previous value until the controlled variable has moved 
slightly beyond the set point. In actual cperation it is the same as hystere- 
sis, aS may be seen in Fig. 2-8. 

A differential gap is caused in the two-position controller of Fig. 2.6. 
If small static friction exists at the bearing on the float arm, the liquid 
level must rise slightly above the desire value to create sufficient buoyant 
force to overcome friction when the level is rising. Also, the liquid level 
must fall slightly below the desired value when the level is falling so that 
the weight force may overcome the friction. This kind of differential gap 
may be caused by unintentional friction and lost motion. 

A differential gap may be intentional, as when a magnet is installed on 
the float arm in Fig. 2.6, causing a hysteresis in float-arm action. Similar 
arrangements are common in domestic thermostats and are employed for the 
purpose of preventing rapid operation of switches and solenoid valves for 
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3. Description of the system, 


The nonlinear feedback control system to be studied is a relay servo 
system as shown in Fig. 3-1. 

Block one contains the transportation lag. The nature of this block 
and its analytic representation will be discussed in section 6. 

Block two is the ideal relay, no dead zone and no hysteresis. 

Block three is the motor load unit or plant. It is assumed to be a 
second order system. In order to have the transfer function always low 


pass in nature, we choose the value of G(s) as 


an oC sae Cap 
Block four is a tachometer. In this system velocity feedback compen- 
sation is used, 
This system, in addition to the velocity feedback, will also use unity 


position feedback. 
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4, Nonlinear analysis. 
4.1 General 

In many cases the system designer is forced to work with components 
which are not entirely linear. In order to linearize the response of a 
nonlinear actuator we use the feedback, The basic difficulty with the non- 
linear equation is that it is not possible to exhibit a general solution. 
Associated with this is the fact that there is no complementary function 
or transient solution made up of terms with characteristic shapes, there 
are, therefore, no system roots Such as were defined for linear systems 
with constant coefficients. It is not the scope of this paper to include 
a comprehensive treatment of nonlinear systems, 

Many syStems are approximately linear when the amplitude of the 
system variables are not too large or too small. For large enough ampli- 
tude all physical systems must become nonlinear. A system which for small 
amplitude is a linear unstable system must become nonlinear when the oscil- 
latory amplitude exceeds some value. Such a system will sometimes reach a 
stable oscillatory motion which is termed a limit cycle. Although such 
motion is usually periodic, it need not be sinusoidal. 

There are other large-amplitude nonlinearities which are in a sense 
destablizing. In these systems, small amplitude response may be stable 
and satisfactory, but large disturbances may send the system into uncontrol- 
led large-amplitude oscillations. 

There are also nonlinearities which become important for small amplitude 
motions, Mechanical friction and backlash are good examples of this type. 
In the case of backlash it is not unusual to find a limit cycle of small 


amplitude which leaves the system capable of following commands or inputs 
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with good accuracy and in an essentially linear mannes except for the small 
limit cycle motion which is superimposed on the desired response. 
4.2 The describing function. 

Many functions involved in automatic control are not desirable by 
elementary techniques, and they must be described by graphical or other 
methods. 

The output of a nonlinear element will not necessarily be sinusoidal 
even though the input is a pure sinusoid, but it will generally contain 
fundamental and harmonic waves. The frequency of the fundamental of the 
output will often be the same as the frequency of the input. [In this case 
a describing function can be used and the element may be treated as a lin- 
ear element, 

Since most elements of a feedback system are of a low-pass nature, the 
higher harmonics may be greatly attenuated in the rest of the loop, so that 
the input to the nonlinear element may be quite accurately sinusoidal. If 
this is the case, it is necessary to compute the transmission around the 
loop for the harmonics. In particular, for the nonlinear element we need 
only a function to describe the fundamental component in the output as a 
function of the amplitude and frequency of the sinusoidal input. This furc- 
tion has been named a describing function, Since the describing functions 
are approximations to the output-input ratio of an element having slight 
nonlinearily, it is anolagous to the transfer function of a linear systen, 
and for certain purposes will be used in a similar manner. There is, of course, 
one important difference. Since a transfer function describes a linear sys- 
tem, it may be frequency-dependent but it is not amplitude-dependent. On 
the other hand, the describing function is usually amplitude-dependent and 


may also be frequency-dependent. The describing function is expressed in 


7, 





terms of the amplitude and phase of the fundamental of the output signal 
related to the amplitude and phase of the input signal. The describing 
function can be obtained analytically or experimentally. 

Describing functions are generally useful in automatic control anal- 
ysis if the neglect of the harmonics in the output signal does not cause 
difficulties. They may be used with adequate results when the element in 
question is followed by one or more elements having appreciable attenua- 
tion of the harmonics in question. 

4.3 Computation of the describing function. 


Consider a nonlinear element whose input is sinusoidal, 
x = X sinat (4.1) 


and whose steady state output can be written as a Fourier series; 


Y = y(t) = a fet a, coscuyt + b, sinawt + a,cos2ut + b, sin 2ut 
+ e#¢6¢6UmG@]hUO 
Ooo 
= an/2 = = (a, coskat + b. Sinko & pee + oe) 
Evidently 
v 
ay = aie f y(t)cos(ka t)d(art) (4.3) 
g 
b = 2/¢ [7 y(t)sin(kes t)d@t) (4.4) 
Oo 


The describing function is defined as the complex number 


@ 


Jf 


N=N, + jN, = NE (4,5) 


1 2 


where 


N, = b,/x N, = a,/X 


or, alternatively, 


i, , wee -j &, 
pies bg ee 


Ne) =GoGw) = OEE — (4.6) 





This definition immediately puts certain restrictions on the permissible 
nonlinearities. The output function must be a periodic function with 
period equal to that of the input. According to this definition, devices 
which produce subharmonic response do not have a describing function. 

The form of Equation 4.1 implies that the input is sinusoidal and 
does not contain d-c bias. If a nonlinearity in a feedback system is not 
symmetric, then in general there will be such a bias. 

Because we are only interested in the nonlinear device of ideal relay 
only ideal relay will be discussed. The characteristic of the ideal relay 


is as shown in Fig. 4-l. 
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Fig. 4-1. Characteristic of ideal relay. 
From Fig. 4-1 
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The describing function of the ideal relay can be obtained very simply 
by carrying the integration indicated in Equations 4.3 and 4.4 and 4.6. And 


in calculating the describing function, we are only interested in the funda- 


mental terms. Thus we have the describing function of the ideal relay as 


Cn es ae 
Ge) = TTX 7 Can) 
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From the above results, we know that the describing function of the 
ideal relay is amplitude-dependent and independent of frequency. This 


function is plotted ite kin Pie. 4.2. 


— 


Fig. 4 Describing function of ideal relay. 

4.4 Analysis of the limit cycle. 
4.4.1 Frequency response method 

Limit cycle angular frequency @, wis predicted by the frequency re- 
sponse method, is such that the loop phase shift is 360 degrees and this 
occurs when the loop frequency response function has no imaginary part. 

The flow diagram for sinusoidal quantities is shown in Fig. 4.3; the 
relay is represented by a describing function in cascade with a time de- 
lay transfer function and the relay output takes on the values + or - unity. 


The condition for 
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a limit cycle is therefore met when 
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by rationalizing this expression the condition becomes: 


SW, Cy 2. : 
e lex, fo, | is real 
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Expanding the exponential term gives the limit cycle condition as 


(Cos ds, t, ~ oe OS, Caan 4 fH) is real 


Equating the imaginary terms to zero yields the condition 
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For a given value of ty (delay time) the dimensionless period P = 
is multivalued. These values correspond to the possible modes of osciila- 
tion and only the longest period is stable. 

According to the frequency response method the servo loop gain at limit 
cycle frequency is unity since the oscillation is exactly regenerated in the 
loop. To determine the limit cycle amplitude at the input to the on-off 
element it is therefore necessary to solve the equation obtained by equating 


the loop frequency response function to unity at the angular frequency, that 


is ; 
“Jol, Cy 


Af zeae =! 


where x is the amplitude of the fundamental component of the error signal 


(4.9) 


when limit cycling. Thus at limit cycle frequency, the imaginary part of 
the loop frequency response function is zero and by rationalizing the left- 


hand side of equation 4.9 and equating the imaginary part to zero we obtain 


Poh ON, Costu ty + Sin @s. td —- 
TX), er 


or 


Xx = <—- a aX Cosa, Cex A See Ld 
“L XG Sep )E cw, ? (4.10) 


where x is the dimensionless measure of misalignment limit cycle amplitude 
and P is the time limit cycle period. 


The values of P and 1 from equations 4.8 and 4,10 are plotted as shown 
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in Fig. 4.4. From Fig. 4.4 we can predict limit cycle amplitude and period 
in an on-off system. 


- 


Equation 4.9 can be rewritten as 


Ce =~ (4.11) 


Equation 4,11 may be analyzed for stability by Nyquies method or by root 


locus method if the describing function G. is amplitude-dependent only. 


D 


When Gy is amplitude- but not frequency-dependent, it is possible to 
reduce a very simple criterion for predicting stability or instability. 
The basis for the criterion is that the leftside of equation 4.11 is only 
frequency-dependent and the right side is only amplitude-dependent. Margin- 
al stability as predicted from the Nyquist criterion and the describing fun- 
ction concept requires a solution of equation 4.11, or in other words, an 
intersection of these two curves. At the point of intersection, the frequency 
parameter on the one curve and the amplitude parameter on the other give the 
frequency and amplitude of the corresponding oscillation. If there is no 
intersection, then the system is presumably stable or unstable for all condi- 
tions of operation. Since it is difficult to find a rigorous definition of 
stability for nonlinear systems, the above statement must be regarded as pro- 


visional, pending more detailed investigation or more careful definition of 


stability for urusual systems, 


eo [kta 
From equation 4,11, we have Cay == Joa(suriy ' 
me TEX 
— ae lj D7es xX 


if we put t, = .2 seconds and by assigning some values of ww , we can have 


d 
the following data; and from these data the Nyquist plot can be plotted as 


shown in Fig. 4.5. 
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Nyquist Plot of Eq. 4-11 
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From the Nyquist plot, a limit cycle is observed; it is at a point where 


Ors 52 and { == 1.71. (The intersection of these two curves). 


intersection may indicate a stability boundary in the sense that excitation 
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at level below that corresponding to the intersection will be unstable for 
the system, and should reach a limit cycle in which the input to the non- 
linear element is roughly sinusoidal with angular frequency G7 =.52 and 
X = 1.71. The system should be unstable for low level of excitation and 
stable for level above that corresponding to the intersection. 
4,4,2 The root locus method 

Because the describing function of an ideal relay is indepeacede 
of frequency, so that the root locus method can be readily applied. The 
nonlinearity introduces no phase shift but only a gain dependent on ampli- 
tude, It will be assumed here that the behavior of the nonlinear system 
can be approximated through use of a root-locus diagram for the system with 


the nonlinearity eliminated, By block diagram manupulation as follows; 


0c . 





Ug Ce) 
Fig. 4.6. 


and algebraic manupulation provides the stability relationship, the 
characteristic equation 


G, G, G,(H +1) = - 1 (4-12) 
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G(S) =- 1 (4.13) 


,, | =-St 
Oe K KS + l)e 
S(Ss + 1 ) 
y, 4y 
Bo 1X 
K. Jel 
For t = .5 second, 
d 
put K. = 1.5 then we have 
: -.58 
G(s) = K (S + 0,67)e | 
7 S (S + 1) 


Equation 4.13 is complex and may be split into two equations by equating 
the magnitude and phase angle on both sides respectively. Thus we have the 


magnitude equation 


[ccs] = 1 (4.14) 
and the phase angle equation 
/G(s) = £ 180° (4.15) 
where 
Lets) = Lots) - ext 
fG(s) -wT = + 180° (4,16) 
where G’ (s) = a 


4 
a method is introduced to obtain the required locus of equation (4.16). 


First it is to construct a family of root loci (or phase-angle loci) of 
the equation 


where §, is a constant phase angle. Thus a family of loci can be plotted 


1 
by various values of 05> Fig. 4.7 illustrates a family of root loci for the 


transfer function G '(s). Since the locus with phase angle 0, is symmetrical 
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to that with phase angle (-@) with respect to the real axis, only the loci 
on the upper half plane are shown, 

Next, it is to construct another family of phase-angle loci (the name 
dead time lag loci is suggested) on the s-plane from the relation 

-WT = O, 

the dead time lag loci family is also plotted in the same sheet as shown in 
Fig. 4.7. Then by superimposing these two families of loci on the same 
s-plane, the point of intersection of the corresponding curves with the sum 
of phase angles 0, and d. equal to 180 degrees are the points of the locus 
required. The curve drawn through these points is the required locus. It 
is the required locus because each point on it satisfies equation 4.16. 
The determination of the gain can be carried out then by the criterion in 
the root-locus method, ‘The magnitude equation is 


1 


= “TG + ja), e 


S 


(4.17) 
which will give the gain K of any point on the required locus. 

The required roots are the points where the value of gain K is equal 
to one. The gain limit, which is the maximum gain below which the system 
can have damped oscillation, can be found readily. This is the value of 
gain at the point of intersection of the required locus and the imaginary 
axis. This value is found to be 1,21. 

This given function transfer has an infinite number of roots, which 
can be obtained by constructing additional branches of required locus. 

In general, the number of influential roots depends on the relative 
magnitudes of the dead time T and the other time constant in G'(s). 

The limit cycle occurs at q@my= 1.67 rad/sec. and k = 1.21 where 


k = k'k 
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5. System with Dead Time. 


5.1. Transfer function associated with pure dead time 
If the input to component with a pure dead time characteristic is f(t) 
u(t), then the output of the component is f(t - T), where T=>O is the 


7 
dead time. The transfer function associated with pure dead time is, there- 


fore, 
-TS 
ie 7 (5.1) 
2 att (5,2) 


For sinusoidal input at frequency a7 , the transfer function given 
in equation 5.2 becomes 
i LJ 
G(jar) = ee I~? (5.3) 
The gain is thus seen to be independent of the frequency and equal to unity, 
and the phase shift is seen to be directly proportional to frequency, 

5.2 Frequency response of a system with dead time. 

It is of interest to compare the frequency characteristics of a unit 


with pure dead time T and a unit with a first-order lag transfer function 


with time constant equal to T, That is, 


G(s) = Se (5.4) 
and 
OG: Se (5.5) 
2 TS +1 : 
The gain associated with G, (jar) and Gi (jar) are 
|S, (iar) = 1 (5.6) 


|G, (ie) = ‘1 + eT Gar) 


and their respective phase angles are 


9. = aT (5.8) 
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and 
0. = - are tangyT (5.9) 
The Bode diagram for the two units are plotted as shown in Fig. 5.1, where 
it can be seen that for wT < 1, the units have nearly the same effect on 
a given input signal. Generally, however, since for higher frequencies the 
dead-time unit exhibits more phase lag and more gain than the first-order 
lag unit, the introduction of the dead-time unit in cascade within a given 
closed-loop system has more adverse an effect on the stability of the system 
than does the introduction of the first-order unit. 

Dead time has the characteristic of a pure phase shift since the magni- 
tude is always one. Seen easily from the graphs drawn in Fig. 5.2 that no 
other dynamic element builds up phase lag so much as dead time. In combina- 


tion with other dynamic elements the phase lag of dead time is simply added 


to obtain the total phase lag. 


Magnitude ratio 
Phase lag 





Frequency, T 


Fig. 5-2. Dead-time element. 
A useful approximation for dead time T is the following function from 


the Pade list, 
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The magnitude of this function is everywhere 1.0 and the phase lag cor- 
responds very closely to dead time up to frequencies given by wT = 0.5. 
In addition, this function may be easier to incorporate into complex system 
analysis. 

The destabilizing effect of the transportation lag, as dead time is 
also known, and clearly illustrated in a Nyquist diagram. For example, if 
a unit with a pure dead-time characteristic given by 

G,(S) = a5 (5.11) 


is connected in cascade with a component with transfer function 


2 
9 = Ser 


as in Fig. 5.3, the Nyquist diagram for the system for various values of 


(5.12) 


the dead-time appears as in Fig. 5.4. As can be seen in this graph, for 
small values of transportation lag the introduction of the dead-time unit 
merely results in a decrease in relative stability, but for large values of 
transportation lag the system can be unstable. 

The effect upon stability of the introduction of a transportation lag 
in cascade within a closed-loop system of the type illustrated in Fig. 5-5 
(a) can be determined from the Nyquist curve for the original system in the 


following manner. The characteristic equation of the original system is 


1 + G,(s)G (s) = 0 (55205) 
and that of the system which is obtained by introducing a transportation 
lag T within either the forward path or the feedback path is 

-Ts 
1 + G,(s)G (sje = 0 (5.14) 


Equation (5.14) can be rewritten as 


G,(s)G (s) eee (5.15) 


or, in terms of frequency response, as 


jatl 


G(54)G, (jar) =e (5.16) 
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Fig. 5-3. Block diagram of a simple feedback 
system with transportation lag. 
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Fig. 5-4, Nyquist diagram for the system 
illustrated in Fig. 5-3. 
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Fig. 5-7 Determination of phase margin in a 
System with transportation lag, 
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The locus of G,(ja7)G, (jar) is the Nyquist curve for the original system; 
equation (5-16) states that the effect of introducing transportation lag T 


is to shift the critical stability point in the GG -plane from -1 + jO to 


b 
-cosa7T - jsinewT. That is, the critical point becomes a critical locus, 


in the form of the unit cycle in the GG -plane, as shown in Fig. 5-6. 


b 
The characteristic equation (5-16) is satisfied only if the locus of 

G. (jar)G, (jar) intersects the stability locus in such a manner that the 
frequency on the Nyquist curve and the frequency on the stability locus are 
the same at the point of intersection. If the two frequencies are not equal 
at the point of intersection, a measure of relative stability is afforded by 
the stability locus. Thus, if the Nyquist diagram for a given system has the 
appearance shown in Fig. 5-7, the phase margin can be determined as follows: 

1. Determine the cross-over frequency GJ. , at whichthe Nyquist 
curve crosses the unit circle. 

2, Locate the point on the stability locus for w=axz. 

3. Draw a radial line segment from the origin through each of 
these points. 

4, Measure the angle between these line segments. This angle is 


the phase margin, with the same positive sense as in the case of a system 


without transportation lag. 
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6. Transportation lag. 


6.1. General 

The phenomenon of dead time, sometimes called transportation lag or 
time delay or ideal time lag, can best be categorized as a lack of change 
in output for a specified interval of time following a change in input, 


this is illustrated in Fig. 6.1. 


e 





Rag, Gl 


Transportation lag may be encountered in process involving pressure, 
temperature, and also may be found in process involving chemical reactions 
when a finite time must elapse before a reaction begins to occur. 

Occasionally in analyzing engineering systems, it becomes necessary to 
Simulate time delay (or transportation lag). For example in studying trans- 
mission processes through hydraulic or pneumatic lines or in studying the 
motion of radio active particles in neuclear reactors it is necessary to 
produce an output voltage which is identical with an input voltage but dis- 
placed in time by a specified amount. 

Transportation lag may be expressed by the following mathematical re- 
lationship: Given f(t) for t =— 0, the function f(t - T) differs from the 
function f(t) by the dead time T, where f(t - T) is defined to be zero for 


eS — T. 
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A direct and accurate method for generating such time delay is to 
employ a magnetic tape recorder in which the space between the recording 
and reading heads is adjustable, by suitably controlling the relative loca- 
tion of the two heads and the speed of the tape a wide range of time delay 
can be simulated. 

Approximate methods which do not require special auxiliary equipment, 
once an analog computer is available, are based upon the Laplace-transforma- 


tion relationship 


Ll fe-W) = Fue = L(fojpe” (6.1) 


-Ss ; 
since e : cannot be expressed as a rational function, the transfer func- 
tion for dead time can be simulated only approximately, thus the remaining 
methods depend on an approximate series or technique of network synthesis. 


The Taylor series expansion of this function is 








-sT a aes) 4 (rst (TS 


( 
ae { — |S + ey 37 Zi! LR! C0) 


The generation of the terms of this series would require successive differ- 
entiation with attendent noise problems, The rate of convergence of equa- 
tion 6.2 is further relatively slow, particularly for large values of T, 
or high frequencies. For this reason the Taylor series is not well suited 
to the generation of transportation lag for problems involving high fre- 
quencies or long time constants. 
6.2 Pade approximation 

A more practical approach to the generation of the transportation lag 
function is to employ the so called Pade’ Bpeeormacions. From equation 


6.1 the transfer function of transportation lag is 
a = (Fe - 7) a _fe (6.3) 
o (f CI] E¢ 
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substituting s = jw to find the frequency response gives 


ae ~Jwy 














——— a = 6 DS mie frre, 
ae 
thus 
22 —/Cos wr + SinraT == / 
EY (6.4) 
= eee gee. yp oe 
DP man ELT ts Kao ae ay) ee aie a (6.5) 


from the above equations the frequency transfer function of the trans- 
portation lag is characterized by an amplitude ratio which is unity for all 
frequencies and a phase angle which varies linearly with frequency over the 
desired range of interest. From equations 6.4 and 6.5 it is seen that the 
quotient of any pair of complex-conjugate polynominals has an amplitude of 
unity and also possess considerable phase shift. 

In the Pade' approximation the exponential is approximated as the ratio 


of two polynominals in s, according to 


ep leaas 7c) a (6.6) 
OP) 


where subscripts m and n refer to highest power of s in the numerator and 
denominator, respectively. Expressed in its basic form, the Pade' approxi- 


mation is 


x E 
eX —— lim Fue (x) (6.7) 
(Usv) —» Co Guy (xX) 
where 
Vex V(v- XO 
oo — | a (CV +v)! + C(UtV)(Urv~1)2 | re 
(6.8) 


ee ViV-Dlv-2)---2 2-/ x” 
(Yt) (Stv-irrres (FI) VY 4 
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(6.9) 





the convergence for this series expansion is quite rapid and often values 


of u and v of 2 will give good accuracy for short time lag and low frequenc- 


ies. 
The Pade’ approximation of eres form = Ve= 2 anduit x = =sl is 
~ ST ,- 267 +367) 
e = 2 i 3 (6.10) 
f/tZ(ST)-+ 3; LS, Ta) 
ees 7 ee Sie 3 
Cc - L (6.11) 


S°T? 4 68ST + 12 

the higher the powers of the numerator and denominator, the more terms of the 
series can be matched and the smaller the error in the time delay function. 
Equation 6,11 is the second order Pade' approximation. The fourth order Pade' 


approximation (u = v = 4) for Sia is 


(ST) BOs T) + 1 §0(S7) - S4.0(ST) +/680 (6.12 
~ (§T)* 4 2067)? + 180 ST)*4+ B40(sT) +1680 aes 
Higher order pade' approximations provide improved accuracy and fre- 
quency response. [It is noted that the order of the characteristic equation 
is raised by the order of the approximant used to simulate ee. If the 
second order approximant is used, two additional roots are introduced into 
the characteristic equation. 


6.3 Simulation of transportation lag 


Equation 5.12 can be rewritten as 











,~ 224 182 _ 8404 1680 

LFC-7T)) ST Sis. GT amor 5) 
APO] ~ = 7, 20 4 {80 , S40 , 1680 
ST (s7= (ST) (8s 7 )F 


since the differential operator, D, and Laplace operator are interchangeable 


A 8 





when initial conditions are zero, equation 6.13 can be expressed as 





i= ES, a ISO _ 840 , 1660 

£(éE-T) ae en Cy ear 
(i 20 , 180 | 40 , 1680 

f is ID T * (p72 (pr)3 * Cam) 





(6.14) 








rearrange equation 6.14 and we get 


fa-1= £(t) - [22 fee) + 22 fe-TIS 
+ Gees na Se Ge 7) | a 
~ [S89 Fee) + FBP £e-7)) Se 
+ [ 4239 Fee) = 1882 4¢-1)) (6.15) 
according to equation 6.15, the analog computer circuit can be simulated 
as shown in Fig. 6.2. 

The time delay T in second is determined by the magnitude of the input 
resistors to the four integrators. Such a circuit works well only ina 
limited frequency range (direct current to approximately 50 cps). In- 
stability results if the variation of the input signal f(t) is too rapid. 
It is therefore not feasible to apply a step function to this circuit. In 
case f(t) is not equal to zero at the commencement of the computer run, it 
is necessary to apply appropriate initial conditions to the integrators. 
The transportation lag devices described thus far are limited to applica- 
tions in which the delay time T remains constant throughout the computer 
run. 

The method of simulating transportation lag by the Pade’ approxina- 
tion yields what may be termed optimally flat phase characteristics at the 
Origin but does not fit the phase characteristics as well over an extended 
frequency range. The major drawback to the use of the Pade' approximation 


to simulate transportation lag is their response to a step input or other 
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input having high frequency components, since the amplitude ratio is 
always unity but the phase shift is approximated only for low qpT pro- 
ducts. The phase shift characteristics of certain of the Pade’ approximation 
are shown in Fig. 6.3. As a result of the inability of this approximation to 
linearly phase shift high frequencies, the output of an approximant rings 
badly for high-frequency inputs. This behavior is illustrated in Fig. 6.4. 
The actual location of the transportation lag transfer function in the 
forward path is immaterial and it may be placed at the end of the path in 
the computer where the high frequency components have been attenuated to the 
greatest possible extent. 
In the great majority of cases the fourth order Pade' approximant is 
more than adequate for the simulation of transportation lag. In cases 
where linear phase shift of higher «aT products than those shown in Fig. 
6.5 is desirable, higher order Pade' approximants or more complex curve 
fitting techniques may be used to simulate yc The phase shift character- 
istics of the sixth order and eighth order Pade’ approximants are shown in 
Fig. 6.5. On the basis of this result it may be inferred that the order 
of the Pade’ approximation may be increased to increase the magnitude of 
the c¢erT product that is linearly phase shifted without distorting the 
phase shift of low “ T products. Because the Pade' approximation may be 
expressed in a general form from which any order approximant may be develop- 
ed, it is felt that the use of curve fitting techniques requires a needless 
expenditure of effort without any tangible benefit. 
6.4 Generation of time delay from magnetic tape recorder. 
The transportation lag generated by using the magnetic tape recorder 
is very accurate, and it saves a lot of time and labor for simulation by 


using amplifiers. 
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The magnetic tape recorder used has the following specifications: 
Ampex Model FR1107 magnetic tape recorder 
Tape speed: 3-3/4, 7-1/2, 15, 30, and 60 inches per second. 
Input impedance: 100,000 ohms unbalanced to ground. 
Input level: 1 volt rms recommended for normal recording level. 
Output impedance: 125 ohms 75 Af, unbalanced to ground. 
Output level: 1 volt rms normal, across 10,000 ohms load impedance at 
normal recording level. 

The response curves obtained from the magnetic tape recorder are as 
shown in Fig. 6.6. Compare this with that obtained from analog computer 
as shown in Fig. 6.4, it is obvious that the magnetic tape recorder is 
most desirable for simulating the time delay, since it gives a more ac- 


curate result. 
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(b) delay time - 1.7 seconds 
Fig.6-6 Curves obtained from magnetic tape recorder 
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7. system simulation. 


7.1. Analog computer simulation 


The block diagram of the system to be investigated is as shown in 


ee 7. le 


b. . 





Fig. 7.1. The block diagram of the relay servo system. 
The simulation of the transportation lag is shown in Fig. 6.2. The simula- 


tion of the ideal relay is as shown in Fig. 7.2. 


Fig. 7.2. The simulation of the ideal relay. 
From the block diagram of the system, the analog computer circuit 
diagram can be simulated as shown in Fig. 7.3. 
Because of its convenience and accuracy, the magnetic tape recorder is 
selected for the generation of transportation lag Ces through the experi- 


ments. 
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7.2. Digital computer program 


The block diagram of the system is shown below. 


Delay 


BCC) tO E(t) 


DA DELAY 





ID2 
ADELAY 


Fig. 7.4. The block diagram of relay servo system. 
The variable names used in the program, E(IT), EDELAY, etc., are shown 
below their mathematical counterparts. Note that the fixed point IT or the 


floating point T indicates the number of increments,starting with IT = 1; 


zero time = DELAT time = 2,DELAT, ete. 
' ae 


| | | | 


i = 1 IT = 2 If =.3 
i= 1. T S20. T = 3. 
E(1) E(2) E (3) 

etc. etc. etc. 


The method of computation is simply to evaluate the analytic expression for 
the response at each increment point. The response is always supposed to 
have started at the time of the last relay operation (TF in increment measure) 
with initial condition equal to the final values of A(IT) = AF, and DACIT)/ 


CTACH = DAF at the end of the previous response segment: 
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new response computation starts here 
| based on relay operation at IT = TF 
lig computation at this point shows change of Sign of EDELAY from previous 


value, the relay operates here, TF, AF, and DAF are recorded as the current 
value of IT, A(IT) and DA(CIT). 





Analytic response. 


The plant transfer function is 


A(s) / 
SiS) St S(s 7%) Cie) 
B(s)== S*A(S) EX Si sy (iy 
and the differential equation describing the plant is 
2 
a ee ae) , ABE) 
t)== “pee +t x “GF (7.3) 


with initial conditions a(0) and a(0) (AF and DAF in program). 


The input is a step V 


4. — 574s) - Sa(o) -d(o) + ASA) - Wa) 


S 
. V-r S7a(o) + S{QAle) +x aco) (7.4) 


i Se) 
This can be expressed into partial fractions 


Bey XZ. XSta)t (star 29 


Se 5 So: an SS Csr oi.) 

- Cie 5) 
a (y+ Z)S t(Xt Ya)S + Xo 
S* (S+2) 
equating the coefficient of s 

S : Yu = a(o) (7.6) 
S: X+Yod== A(0) + dao) (747) 
os nn (7.8) 
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kK == 

ae 
Y= 2] Q(o)+ dae) - X} 
Z=ace)- 7 


from equation 7.5, the time solution is 


V/ALPHA in program 


(DAF + ALPHA-X)/ALPHA 


AF - Y 


il 


AO= NG aa +ze" 
= X*TAU + Y + Z*EXPF(-ALPHA*TAU) in program 
A(t)=X-szae** 
= X - Z*ALPHA*EXPF(-ALPAH*TAU) in program 
t-~-—> TAU is time from last relay operation 
Program accuracy. 

1). The time delays can be introduced only as an integral number 
of increments (each of magnitude DELTAT). 

2). The only actual inaccuracy is the fact that the relay can only 
operate at integral values of IT*. The accuracy, for a given increment size, 
DELTAT could be fairly easily improved by introducing an interpolation pro- 
cedure between computation points. The same technique could be used to 
introduce non-integral number of the delays. 

The plant dynamics can be easily changed by re-working the partial 
fraction expression and the expression for a(t) and a(t). 

The program for the relay servo system is shown below. And one case 


of its output (t, = .8 second, k = 2.0) is attached to this program. 


*This error may be regarded as a variability in the forward path time delay, 
so that the actual delay may lie between the desired ID1l and ID1+1l. 
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7.3 Comparison of the results obtained from the digital computer and from 
the analog computer. 
The result obtained from the digital computer for the system with time 
delay .4 second has been checked with analog computer. They are quite in 
agreement with each other as shown in Fig. 7.5 and 7.6. 


comparison of the results 


by digital by analog 

computer computer 
overshoot and undershoot 35 2.6 
rise time 5.9 5.7 


Because the digital computer can give an accurate and exact value for 
any desired instant, therefore, the digital computer is selected for con- 
ducting most of the experiment work, and the analog computer is used to 


check the results obtained by digital computer, They are quite in agreement. 
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Fig. 7-5 Output,velocity of the systen, ts - -4 second 
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Fig J — ; Phase trajectory of the system with ts - -4 second 





8, Effect of transportation lag on relay servo system. 
8.1. General 

From the proceeding sections, we know that the time delay has an adverse 
effect on control systems, and which is unavoidable in nearly all the control 
Systems. Because the time delay may basically be represented by an infinite 
series of terms, this creates an infinite number of roots. 

Later on in this section, it will be noted that the time lag results in 
successive overshooting and undershooting of the correct final value, re- 
presenting serious oscillations about the final position. In addition to 
this it will also change the root locus and phase trajectories, 

8.2 The frequency response. 

From Fig. 4.6(c) we have the transfer function with Gh removed 

= Sas 


a Clio = 
GISJ== —“Stse1) (8.1) 


where KR and T are variables. In case of time delay equal to .5 second, it 
is found that the tachometer gain, ky = 1,5, will compensate the system to 
an acceptable optimum operation, this will be seen'’in section 9, now equation 


8.1 can be written as 


~~ 


.. Vara oe 
Gb)—= rr ae (8.2) 


The Bode diagram for the above transfer function is shown in Fig. 8.1. So 
far as the linear system is concerned, the gain margin foe »2 second time 
delay system is approximately 1.6 db, but for the undelayed system, it is 
stable for all gains and never unstable. 

The time delay has introduced additional branches into the diagram. 
These branches appear due to the periodic nature of the time delay, and the 


increasing phase shift introduced at higher frequencies. Yt is evident that 


72 





"7-8 °by FO weaseltp spog [-g °31y 











| 
2 | 
| | 

t 

$ 


09- 


aVE€- 


\ 


Do 





the gain margin for each of the branches increases. This is analogous to 
the statement that the roots of the primary branch would go unstable first 
if the gain were raised. 

The gain vs. logur curve of the delayed system is identical to that 
of the undelayed system, since the time delay introduced no additional gain 
to the system. The phase angle vs. logger curve of the delayed system is 


drastically changed since the equation for the phase angle is 
Pp mee O- ae CCe a) ie oe 


The eT term in the above equation has a predominate effect to the phase 
angle, it introduces additional phase shift and extra branches. 

The limit cycle occurs at cy =1.7, k = 1.20. 
8.3 The root locus. 

The root locus of any system will be greatly changed when a transporta- 
tion lag was introduced. As a matter of fact, for an undelayed system the 
transfer function in Equation{(8.2) becomes 


15 (S +°67) 


Seam (8.4) 


no limit cycle exists, because the root locus never intersects the imaginary 


axis as shown in Fig. 8.2, which shows a stable system for all gains. 


Je 





Fig. 8.2.  Reot lieeus of Equre.4. 


But when time delay is introduced this situation will be changed. The root 
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locus for Eq. 8.2 (.5 second delayed system) is shown in Fig. 8.3. The 
construction details for this root locus were described in section 4.4.2. 
Fig. 8.3 illustrates the effect of time delay on the root locus. It 
is apparent that this system will go unstable with only a smali increase 
in gain. The primary branch stability limit before the root locus corsses 
the imaginary axis are k = 1.2., and « = 1.76 rad/sec., where the limit 
cycle occurs, this almost checks with Bode diagram where k = 1.20, aw=1./7 
rad/sec. 
Because time delay introduces periodically phase shift of 360 degrees, 
so there presents additional (upper) branches when at high frequencies. It 
is noted that for the same gain the additional roots on the upper branches 
are attenuated more than the roots on the primary branch. This substantiates 
the Bode diagram findings, which shows that the roots on the primary branch 
phase shift would cause instability. 
8.4 Transient response. 
Because of time delay introduced in the system, there is no graphical 
method that can be easily used to produce the transient response of any 
system. The transient response of the delayed relay servo system (second 
order) for various delay time and various ke (tachometer: gain) are shown 
in Fig. 8.4. In Fig. 8.4(a) there are six curves under the same tacho- | 
meter gain (k, = .02); curve (1) is an undelayed system, the rest of the 
five curves are with .1, .25, .5, .8 and 1.0 second delay time respectively. 
For 1 second time delay system it presents an overshoct and undershoct of 
approximately 60 percent. 
Table 8.1 is the data obtained for ste : i t under different values 


t | 
of ke and ta Plotting these data into curves it is found that they appear 
as a straight line. In Fig. 8.5 are shown and a versus delay time; in 

C 
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Fig. 8-3 Root locus of relay servo system with time 


delay of .5 second. 
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Table 8,1. Data for first overshoot 7 and i. under various time delays 
7 t 
for a step input. 


K. Ty os se amplitude of 
t limit cycle 

02 20 i285 3.50 2.828 no limit cycle 
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: 2a e250 4,10 3,095 1625 

i a, |e 72510 4,60 35.390 2579 

H 8 L500 5. 10 3.650 4910 

A Us 1.6500 5.50 3.850 » 5900 

ae 50 b.0s70 32.05 2.900 no limit cycle 
4 mal 1.0870 3251 2.950 negligible 

B25 1.1630 Sele S060 057 
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+ Wo 1.0550 4,10 7 0 2137 
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Fig. 8.6 is shown the rise time versus delay time, there it shows that the 
effect of time delay on the transient response is to increase the over- 
shoot and undershoot and cause continuously oscillations and make the rise 
time longer, it turns out that the response of the system can never be 


settled. 

It is noted that M ; oe t increase as the delay time increases, 
But when we increase POR of ke they decrease, And to a certain value 
of ke some of them become acceptable optimum operation this will be seen in 
section 9.2, 

In Fig. 8.7 are shown the plant output, rate output and system error. 
Fig. 8.7(a) is the response of the undelayed system, the plant output set- 
tled at 7 seconds; Fig. 8.7{b) shows the plant output, rate output and 
system error for a system with delay time of .5 second and ky, = .02,/ tne 
output continuously oscillates with an overshoot and undershoot of approxi- 
mately 30 percent; Fig. 8.7(c) shows the plant output rate output and system 
error for a system with time delay of one second and ke = .02. The plant 
output continuously oscillates with an overshoot and undershoot of approxi- 
mately 60 percent, the frequency of oscillation is slower than that of .5 
second time delayed system; Fig. 8.7(d) shows the plant output, rate out- 
put and system error for a system with delay time of 2,5 seconds and ke = 
.02, the plant output continuously oscillates with an overshoot and under- 
shoot of approximately 135 percent, the frequency of oscillations is much 
slower than that of .5 second time delayed system, 

8.5 Time delay in the phase plane diagram. 

The effect of time delay in application or removal of the applied 
torque after the control has signaled that the operation should be perform- 
ed. It will be of interest to study the effect in terms of the phase plane 
diagram. 
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Because of the time delay the trajectory in the phase space will be 
prolonged beyond the point c’ = R’ by an amount such that the difference in 
coordinates is just equal to the time delay. That is, in Fig. 8.8 the torque 
should reverse at a, but it does not reverse until 6B is reached, and (Cc. ~ CA 
(ve = v,/2 is equal to ty » Similarly, the torque is reversed at point D 
instead of at c, and at F instead of at E. [t is apparent that a curve could 
be drawn through points such as B and F so that the switching of torque 
direction would take place when the trajectory crossed this line, and the 


line from B to D would play the same role as the line through R' for the 


control without time delay. 





Fig. 8.8. Normalized phase plane diagram for relay servo 
system with a time delay of .l second, for a 
Step unit. 
The effect of time delay is thus to skew the switching line from the 
verticle in the phase space, and quite approximately the angle of skew is 
-tan 't, and is in clockwise direction. 


It is seen in Fig. 8.8 that the trajectory appears to repeat itself, 


since the last two oscillations are nearly identical. The control oscillates 
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with a constant anplitude on the trajectory, Suth a repetive diagram re- 
presenting the steady state oscillation is called Limit cycle, From the 
experimental work it shows that the relay contrcls with time delay and ne 
dead zone exhibits the defect of continuous oscillations. 

tn Fig. 8.9 is shown a family of trajecturies so arranged that the 
initial and final points of each are symmetrically located with respect to 
the origin. Thus each trajectory represents just such 4 limit cycle, and 
Since the time delay for each is fixed, the lecus of the end points of these 
limit cycles measures the time delay required to »sbtain the amplitude of 
oscillation. As would be expected, the larger the time delay the larger is 
the magnitude of the oscillation, and with that, the lower is the frequency 
of oscillation as indicated in Fig. 8.5. In Fig. 8.16 is shown a plot of 
time delay and oscillation amplitude. It turns out to be a straight line, 
therefore, we can obtain an equation to express the relation of oscillation 
amplitude with time delay as 


A OSC 


the above equation is for a second order system specified by 


| aoe / 
G(s) = BSCS ee 


It is interesting to see the relay type control that can continuosuly 





eae one (8.4) 


cycle when a time delay is present. The relay type of control always 
oscillates about the equilibrium position. For any time delay the system 
would always cscillate, and with a magnitude of oscillation depending on 
the time delay. The relay type control can never be umstable in the sense 
that the amplitude of oscillation will continue to increase without limit. 
It is obvious that this is the advantage of the relay servo system, 

In Fig. 8.11(€a) the phase plane of the relay servo system with a time 


delay of zero and .5 seconds is plotted on the same sheet. In Fig. 8.11¢b) 
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is the phase plane of the relay servo system a time delay of one and 2,5 
seconds plotted on the same sheet. From these phase portraits we find a 
most important aspect of the phase portrait is the switching line which is 
not vertical. The effect of the dead time is to rotate the switching line 
into the first and third quadrants, the result is zrowing oscillation. 

The process of moving the switching line from vertical axis to a line 
in the first and third quadrants is seen to affect the nature of the response 
very unfavorably. On the other hand, if the switching line can be made to 
lie in the second and fourth quadrants by means of phase advance, or lead, 
on the error signal, the relay servo system is likely to be a very well be- 
haved one, in that it has small overshoot Pie will reach and remain within 
a given value of error in a short time. This will be seen in section 9. 

Dead time will result in a terminal limit cycle which may or may not be 
acceptably small, the amplitude of oscillation depends upon celay time as 
shown in Fig. 8.9 and .10. 

8.6. The effect of time delay in steady state oscillation frequency. 

In Fig. 8.12 the curve of system error against time is shown. [nitial 
condition of the error is due to a suddenly applied step input signal, this 
gives an initial error and velocity (output rate) to the form of response 
waveform as shown in Fig. 8.13. After the transient has died out the out- 
put shaft follows the input by performing oscillation about a final posi- 
tion as shown in Fig. 8.5. 

In Fig. 8.12 curves 1, 2, 3, 4, and 5 are for a delay time of .1, .25, 
.9, .8 and 1.0 second respectively. [t indicates that the smaller time 
delay causes high frequency oscillation with smaller amplitude, in other 
words, the frequency of the oscillation is higher, and the larger value 


of time delay results in larger cscillation amplitude with lower frequency 
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as shown in Fig. 8.12. 
In Pie. 6.13 cimveso 12. 2.94, 5, “andubGsatettco: a delay time of 0, .1, 
w2D, .2, +8, and 1.0 second respectively. Th® curves show the relation of 
the rate oufput to time delay. 
From Fig. 8.5 or Fig. 8.12, the data for oscillacion frequency against 


time delay can be obtained as shown in Table 8.2. 


Table 8.2. Nata for steady state oscillation frequency vs. time delay. 


delay time Frequency 
ok 2455 
e205 294 
re, - 200 
28 . 156 
iL8) S130 


The data plotted in log log scale sheet there appears to be a straight 
line as shown in Fig. 8.14; it is seen that the frequency decreases as 
the delay time increases. Because the curve of oscillation frequency 
against time delay is a straight line, therefore we can obtain an equation 


to express their relation as 


i oof Diy 
f ni eo (8.5) 


From Fig. 8.12 we take the steady state error for each time delay, 


tabulate the data as shown in Table 8.3. 


Table 8.3. Steady state error vs. time delay. 


Delay time Error 
0 0 
al 709 
we 218 
3 oe 
aS. 484 
wae » 388 


Fig. 6.15 is the plot of error vs, delay time, [t comes out approximately 


~ 


a straight line, it indicates that the steady state error is proportional 
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Figs 8-13. Velocity vs.time curves showing the effect of time delay on velocity. 
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to the time delay, 
§.7. The effect of the magnitude of the relay voltage. 

From the experimental work, we find that the magnitude of the relay 
volteze also has sowe effects on the response curve and phase trajectory of 
the systea., When the relay voltage is smaller the response has a longer 
rise time and a smaller overshoot and undershoot, and there is presented 
some velocity saturation in the phase trajectory, the amplitude of the 
limit cycle is smaller. 

In Fig. 8.16 shows two curves with the same delay time of .5 second 
and same value of Kk = ,02, but with different relay voltage, curve (1) is 
with V =e2 volt and curve (2) is with V =«¢5 volt. From these curves it is 
indicated that the smaller relay voltage causes smaller oscillation amplitude 
and lower frequency, but it presents a longer time rise which is its dis- 
advantage; so far as the time rise is concerned the larger value of relay 
voltage is desirable, therefore, in all the experiments we set the relay 
voltage ecual to .5 volts. The rise time for .2 volts of relay output is 
6.5 seconds, and the rise time for .5 volts is 3.35 seconds. The overshoot 
and undershoot for .2 volt system is 12,5 percent, and that for .5 volt is 
30 percent. 

in Fig. 8.17 the phase plane with the same delay time and same value of 
ke is shown, but with different relay voltage, it presents a velocity satura- 
tion for the system with .2 volts of relay output, and the limit cycle is 


smaller, 
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9. Elimination of the steady state oscillation 


In many applications steady state oscillations of the output shaft are 
inadmissible; besides causing continual operation of relay contacts and 
wear of gears and bearings the power consumed in oscillations may be objec- 
tionable. The simplest method of eliminating the oscillation is to intro- 
duce a dead zone in the relays so that the error always returns to the 
dead zone after a disturbance. But the relay used has been thus far assumed 
an ideal relay with no dead zone, so we shall not discuss this possibility. 
9.1 Velocity feedback stabilization 

A very effective method of reducing the settling time and limit cycle 
oscillation of an relay-servo is to make the relay switching conditional 
not only on misalighment but also on its rate change. This can be achieved 
by feedback within the position control loop if the feedback signals can 
be added algebraically to the misalignment signal at the error sensing 
device. 

The effect of negative velocity feedback is illustrated on the phase 
plane diagram in Fig. 9.1; and the corresponding time response is shown 
mmerig., 9,2, 

It is obvious that the larger the value of k. the smaller the ampli- 
tude of the limit cycle, and the more the switching line shifted to counter- 
clockwise direction; and the amplitude of oscillation is reduced; the over- 
shoot and undershoot for ke = .02 is about 15 percent, for kK =o. 1S 
about 2 percent as shown in Fig. 9.2. 

9.2 Optimising the transient response. 
Because of the introduction of the time delay in the system, it is 


difficult to derive a method to determine the value of ke which will give 
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an optimum operation, which can only be determined through experimentation, 
As we know the system damping will be improved by the tachometer feed- 
back, the net effect of the derivative signal is to cause the relay to 
operate early, so by using kK. we can control the switching characteristics 
of the relay. For many applications of relay servomechanism, the tacho- 
meter feedback compeneation is adequate. 

Through experimentation it was found that the system with a delay time 
below .5 second can be compensated for 5 percent overshoot and undershoot 
which is considered to be acceptable. 

The response and phase plane of the compensated system for various time 
delays are shown in Fig. 9.3(a) through mee, $.3{e),. They have been compen- 
sated by velocity feedback, The values of ky (which give acceptable optimum 
operation) are as follows: 


Table 9.1 Optimum value of k, for eliminating the effect 
of time delay. 


t k o/o overshoot 
and undershoot 


0 23 / 
ul »4 a 
eo 8 1.6 
2 aD ee 


For systems with .8 and 1.0 second time delay, it is found that they can 
not be compensated for 5 percent overshoot and undershoot, even when we 
increase the tachometer gain, Key to a considerable value, the amplitude 
of steady state oscillation is still rather large as shown in Fig. 9.4 and 
Fig. 9.5; the only effect of increasing ke at this time is to shift the 
whole response curve downward, as shown in Fig. 9.5 the response of the 


System with one second time delay at ke = 4,5 is rather low but still with 
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considerable oscillation amplitude which is almost the same as that when 
kK = 1.5, which is beyond the acceptable value. 

Fig. 9.6 illustrates the time delay vs. k, curve at the range from 0 to 
.9 second time delay. The curve appears to be a straight line, but when 
the time delay above .5 second occurs this situation does not hold, since 
the system with a delay time above .5 second cannot be compensated. When 
the time delay is less than .5 seconds, the system can be compensated; and 
because t, vs. k. curve is a straight line, therefore, it can be expressed 


d E 


by an equation as given in equation 9.1. 
Kya 2 + 28, (9.1) 


Table 9.2 and 9,3 are the results for comparing the uncompensated and 
compensated system, The data in Table 9.2 is for the undelayed system 
for which k. = .3 will compensate the system to give optimum operation. 
The data in Table 9.3 is for a .5 second time delayed system for which k. 


= 1,5 will compensate the system to give acceptable optimum operation. 
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iC, The response of relay servo system to ramp function input. 


In the previous sections we have been concerned primarily with the 
response to step function input. Experimental results obtained from the 
relay servo system with a ramp function input are shown in several graphs. 

When a unit ramp input is applied to the relay servo system, the out- 
rut ane the error of the system are shown as a straight line, and the velocity 
for most part is a horizontal line, as shown in Fig. 10.1. The Cdot vs. © 
and Edct vs. E phase planes have the same shape, merely the Cdot vs. C plane 
appears above X- axis. Fig. 10.2 shows the Edot vs. E plane. This is with 
a relay voltage equal to .5 volt. When the relay voltage is increased to 
1.0 volt the same situation remains. For relay voltage at .5 and 1.0 voit 
the only effect of the time delay on the relay servo system is to cause some 
delay in the response, no oscillation has been observed. 

But when the relay voltage is increased to 1.5 volts, this situation is 
changed. For an undelayed system the output still re as a straight line 
but the system error is no more a straight line, and the velocity curve is 
with some oscillations as shown in Fig. 10.3. Cdot vs. C and Edot vs. E 
phase planes have a different shape as shown in Fig. 10.4 and Fig. 10.5. 

For delayed system the output no more remains a straight line but with some 
oscillations, and the system error and velocity are much different than that 
of the system with relay voltage of .5 or 1.0 volt as shown in Fig. 10.6. 
The Cdot vs. C and Edot vs. E phase planes are quite different as shown in 
Fig. 10.7 and Fig. 10.8 respectively. Oscillations and limit cycle has been 
observed when the relay voltage is 1.5 volts. 

The effect of time delay on the response of the system can also be com- 
pensated by velocity feedback. For the .5 second time delayed system, when 
k, = 2.0, the response can be compensated to an acceptable optimum operation 


c 


as shown in Figures 10.9,16, 11. 
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Fig. 10-2, Phase trajectory of relay servo system, relay voltage 
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Ld, dGone lusiions. 

The aim of this thesis was to demonstrate the effects of transportation 
lag on the relay servo system, 

For conducting the experiment by means of an analog computer, first 
could simulate the time delay by Pade's approximation or by magnetic tape 
recorder. [It was seen that the generation of time delay from magnetic tape 
recorder presents a more accurate and better result than that from the Pade’s 
approximation. Therefore, the magnetic tape recorder was selected to 
generate the time delay throughout the experiment. 

In this investigation most of the experiments were conducted by using 
a digital computer which gave more accurate results than the analog com- 
puter does; especially since it can give an exact value at any desired in- 
stant. The analog computer has also been used to check the results ob- 
tained from the digital computer. It is noticed that they are quite in 
agreement. 

From the analysis and experiments the following conclusions has been 
obtained. 

The root locus of the system was greatly changed due to the intro- 
duction of the transportation lag; and because time delay introduced 
periodicaly a phase shift of 360 degrees, it results in an additional (upper) 
branches when at high frequency. 

The effect of time delay on the transient response is to increase the 
overshoot and undershoot and causes oscillations and also leads to a longer 
rise time, the response of the system can never be settled. 

The effect of time delay on the phase plane is to rotate the switching 


line into the first and third quadrants, the results is growing oscillations 
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In other words the transportation lag will result in a terminal limit 
cycle which may or may not be acceptabiy small. 

The plot of time delay versus oscillation amplitude came out to be a 
straight line; the delay time multiplied by a constant of .6 gives the 
magnitude of the oscillation amplitude. 

The time delay versus steady state oscillation frequency when plotted 
on leg log coordinates came out to be a straight line in such a relation 


that a constant of .135 divided by ee gives the oscillation frequency. 

The steady state error and the velocity are proportional to the eine) 
delay. 

These effects caused by the transportation lag can be compensated to 
some extent. The method used in this ees ciracien to eliminate steady 
state oscillation and to optimise the transient response by employing tacho- 
meter feedback, It was seen that a system with a delay time below .5 second 
can be compensated to an acceptable small oscillation amplitude (5 percent). 
But for time delay above .8 second it cannot be compensated to such a small 
amount; the only effect of increasing tachometer gain is to shift the whole 
response curve downward with the same amount of steady state cver- and under- 
shoot. 

The system with a time delay below .5 second can be compensated to reach 
an acceptable optimum operation. For delay time below .5 second its relation 
with tachometer gain is approximately a straight line, therefore, it can be 
expressed by an equation. 

When a time delay is present the relay type control can continuously 
oscillate, and with a magnitude of oscillation depending on the time delay. 
The relay type control can never be unstable in the sense that the amplitude 
of oscillation will continue to increase without limit. This is one advén- 


tage of the relay servo system. 
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